Neel order, ring exchange and charge fluctuations in the half-fifled Hubbard model 



o 
o 

(N 



X3 
O 

o 



> 

cn ■ 
cn : 
o . 

o 



o 



X 



J.-Y. P. Delannoy/'^ M. J. P. Gingras,^ P. C. W. Holdsworth,i and A.-M. S. Tremblay^ 

^ Laboratoire de Physique, Ecole normale superieure de Lyon, 46 Allee d'ltalie, 69364 Lyon cedex 07, France. 
''Department of Physics, University of Waterloo, Ontario, N2L 3G1, Canada 
Departement de Physique and RQMP, Universite de Sherbrooke, Sherbrooke, Quebec, JIK 2R1, Canada. 

(Dated: February 2, 2008) 

We investigate the ground state properties of the two dimensional half-fiUed one band Hubbard 
model in the strong (large- 17) to intermediate couphng limit (i.e. away from the strict Heisen- 
berg limit) using an effective spin-only low-energy theory that includes nearest-neighbor exchange, 
ring exchange, and all other spin interactions to order t{t/U)'^. We show that the operator for 
the staggered magnetization, transformed for use in the effective theory, differs from that for the 
order parameter of the spin model by a renormalization factor accounting for the increased charge 
fluctuations as t/U is increased from the t/U — » Ifeisenberg limit. These charge fluctuations lead 
to an increase of the quantum fluctuations over and above those for an S = 1/2 antiferromagnet. 
The renormalization factor ensures that the zero temperature staggered moment for the Hubbard 
model is a monotonously decreasing function of t/U, despite the fact that the moment of the spin 
Hamiltonien, which depends on transverse spin fluctuations only, in an increasing function of t/U. 
We also comment on quantitative aspects of the t/U and 1/S expansions. 



I. INTRODUCTION 

Effective low-energy theories are constructed and used 
in essentially all fields of Physics. The exponential re- 
duction in the size of the Hilbert space that occurs in 
such theories generally offsets the disadvantage of work- 
ing with the non-local operators induced by elimination 
of the high-energy states. In the context of strongly- 
correlated electrons, spin-only Hamiltonians, such as the 
Heisenberg model, are examples of effective low-energy 
theories that apply when interactions are very strong. 
Present research interests focus on systems in the strong 
to intermediate coupling regime, where one might ex- 
pect that weaker interactions lead to increased electron 
mobility, which in turn should reduce the stability of 
magnetic phases. In the effective Hamiltonian, the in- 
creased electron mobility is taken into account perturba- 
tively by including increasingly higher order corrections 
to the effective low-energy theoryii?''^. More specifically, 
the effective low-energy spin Hamiltonian, Hs, derived 
from the Hubbard model, away from the strict Heisen- 
berg limit [t/U 0], contains conventional Heisenberg 
pairwise spin exchange as well as so-called ring (or cyclic) 
exchange terms that involve n— spin (n > 2) interac- 
tiongSil. These corrections alter the low-energy excita- 
tions and theoretically they may, if large enough, pro- 
duce exotic ground states^. At the present time, there 
are still many experiments designed to search for evi- 
dence of ring exchange terms in materials such as parent 
high-temperature superconductors^ . 

There exist several methods for deriving effective low- 
energy theories such as degenerate perturbation theory, 
canonical transformation, resolvent and projection ap- 
proaches. Their equivalence has been recently demon- 
strated"^. Among the most widely used is the so-called 
canonical transformation (CT) based on original ideas of 
van Vleclii. The main idea behind the CT is the fol- 
lowing. In performing a CT, the "excursions" of the 



degrees of freedom outside the low energy subspace are 
taken into account in the effective theory by non-local 
effective interactions. The true ground state eigenvec- 
tor of the theory is in essence "rotated" to lie in the 
ground state of the subspace of the effective low en- 
ergy theory. The CT method has been extensively used, 
even outside condensed matter physics. For example, 
Foldy and Wouthuysen employed the CT to derive the 
lowest order relativistic corrections to the Schrodinger 
equation, starting from the Dirac equation^. Two of the 
best known applications of the CT in condensed matter 
physics are the derivation of the Kondo model from the 
Anderson impurity model using the Schricffer- Wolff CT— 
and the derivation, starting from the Hubbard model of 
an exchange spin Hamiltonian with ring/cyclic exchange 
termsi^. One important technical aspect arising in the 
construction of effective theories is that operators defined 
in the bare high-energy theory must also be canonically 
transformed before they can be employed in calculations 
within the effective low-energy theory. In the context 
of condensed matter systems, the importance of prop- 
erly transforming operators in high order perturbative 
approaches, used to eliminate states from the high energy 
sector of the theory, has been emphasized in a number of 



situations 
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Consider the Hubbard model with the two energy 
scales defined by t and U, where t is the nearest-neighbor 
hopping constant and U is the on-site Coulomb en- 
ergy. In an effective low-energy theory derived from 
the Hubbard model and limited to the spin-only sub- 
space, the electron hopping processes beyond nearest- 
neighbor lead to a 4-spin ring exchange term, Jc and 
to second and third neighbor exchange interactions, J2 
and J3, which are all of order {t/UY smaller than the 
nearest- neighbor exchange Ji « AP/U. Several recent 
studies have investigated the effect of Jc on the proper- 
ties of a 5 = 1/2 nearest- neighbor Heisenberg antifer- 
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In two dimensions (2D) it is found 
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that introducing a small Jc initially decreases the quan- 
tum fluctuations of the Neel order parameterii^ti^. Simi- 
larly, in a one-dimensional {ID) two-leg ladder the spin 
gap decreasesi2*i^ and the staggered spin— spin correla- 
tions increasei^ as Jc is first increased, again indicat- 
ing a reduction of quantum fluctuations. These studies 
consider as a phenomenological parameter in a spin 
model without reference to the microscopic origin of Jc 
from a Hubbard-like model. However, a tempting inter- 
pretation of the above results is that an increase oit/U 
away from the Heisenberg limit increases the Neel order 
parameter, M^. This picture is re-enforced by a recent 
self-consistent Dyson-Maleev spin- wave calculationi^ us- 
ing Hs derived from the Hubbard model to order t{t/U)^ 
(see Refs. TT^). It is found that Aft for < t/t/ < 1 is 
increased above the value for the Heisenberg limit^^. This 
ensemble of results for effective theories suggests that AI^ 
should pass through a maximum value at some finite t/U 
— a conclusion which is difficult to understand on phys- 
ical grounds and inconsistent with Hartree-Fock calcula- 
tions^*' and series expansionSi, as well as exact diagonal- 
ization of the Hubbard model on small clusters^^ where 
the structure factor measuring staggered magnetic cor- 
relations monotonously decreases as t/U increases from 
the Heisenberg limit. 

In this paper, we show that the above paradoxical in- 
crease of as t/U increases in spin-only theories when 
compared to calculations on Hubbard modelsSfliSiiS^ is 
simply resolved if the magnetization operator is trans- 
formed using the unitary transformation, that eliminates 
order by order the double occupancy states in the Hub- 
bard model^'?*^*i2iii. As t/U is increased beyond the 
Heisenberg limit, so that ring-exchange terms have to be 
taken into account, one can identify two quantum correc- 
tions to the antiferromagnetic Neel order parameter. The 
first correction is an overal amplitude renormalization 
factor, coming from short-range charge fluctuations. The 
second is the usual transverse spin fluctuations. The am- 
plitude renormalization factor from charge fluctuations 
can be obtained exactly and does not need further calcu- 
lation, while the spin fluctuations can be taken into ac- 
count to a good degree of approximation using the most 
naive application of the usual methods, such as the \/S 
Holstein-Primakoff method. The i/C/ dependence of the 
two effects go in opposite directions but the amplitude 
renormalization factor coming from charge fluctuations 
dominates, resolving the above paradox and giving the 
physically correct trend of a decreasing Neel order pa- 
rameter as t/U increases from the Heisenberg limit. An 
analogous result was found in Ref. where it was ob- 
served that finite t/U corrections lead to an increasing 
Drude weight as t/U increases, in contrast to the t — J 
model, where the weight decreases as t/U increases. 

The paper is organized as follows. Starting from the 
Hubbard Hamiltonian, we recall in Section IIA the key 
steps that lead to an expression of an effective spin Hamil- 
tonian, based on the CT methodii^. We then apply this 
method in Section HE to derive an expression for the 



magnetization operator M''' in the low energy, spin only, 
effective theory. In order to expose the quantitative im- 
portance of the charge-fluctuation-induced renormaliza- 
tion factor acting on M\ we present in Section HI re- 
sults from exact diagonalization on small clusters, and 
l/S spin wave calculations in the thermodynamic limit. 
We end the paper with a brief conclusion in Section IV. 
An appendix gives some of the details used in passing 
from the fermionic description of the effective theory to 
the SU(2) spin only description. 



II. STAGGERED MAGNETIZATION 
OPERATOR IN THE LOW-ENERGY THEORY 

A. Effective Hamiltonian 

We begin with the one band Hubbard Hamiltonian, 



uY^n. 



(1) 



The first term is the kinetic energy term that destroys 
an electron of spin a at site i and creates it on nearest- 
neighbor site j. The second term is the on-site Coulomb 
interaction : it costs an energy U for two electrons with 
opposite spins to remain on the same site i] and where 
= cJ^Ci^CT is the number operator at site i. We de- 
rive the low-energy theory using the canonical transfor- 
mation method first used by Harris et al} in this context 
and applied to higher order by MacDonald et ali. The 
method introduces a unitary transformation, e'"^, that 
"rotates" ifn into an effective spin-only Hamiltonian, ffg, 
and corresponding state vectors into the restricted spin- 
only (SO) subspace. 

As introduced in Refs. |ll2j |. the transformation relies 
on the separation of the kinetic part T into three terms 
that respectively increase by one (Ti), keep (Tq) or de- 
crease by one (r_i) the number of doubly occupied sites. 
We write: 



= Ti 



To + r_ 



where 



(2) 



(3) 



Tn = -t 



(5) 



where a stands for up if a is down and for down if cr is up. 
This decomposition of T comes from right multiplying 
the kinetic term in Eq. ^ by rii^^ + hi,g — 1 and left 
multiplying by nj.g. -|- hj^g — 1. 
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Using the HausdorfF formula, the transformation e'*^, 
appUed order by order in t/U to i?Hj gives: 



n - .^Srr _ rr , I*'^, Ha] , [iS , [iS , Ha]] 

J^s — e Hue — Mn H 



1! 



2! 



This unitary transformation is chosen so that, to a given 
order, the resulting Hs does not change the number of 
doubly-occupied sites. Order by order the weight of the 
states with double occupancy are reduced. For the com- 
plete transformation, Hg and the corresponding ground 
eigenstate vector |0)s are completely confined to the SO 
subspace. 

To third order in t/U and following MacDonald et al? 
we recover the expression of the generator iS^^^ of the 
unitary transformation e*'^' ' as being: 



1 

U 
1 

IP 
1 



1 

CP 



([Ti,To]-[ro,T_i]) 



(- [To, [To, m - [To, [Ti, To]] - [Ti, [Ti, To]] 
2 



[r_i,[To,r_i] 



-[r_i,[Ti,r_i] 



[Ti,[Ti,r_i] 



(7) 



This expression combined with ^ leads to the expression 
of the third order expansion of the Hubbard Hamiltonian 
in terms of the operators Tm introduced above. Defin- 
ingS T'^'')(mi,TO2, • • • ,TOfc) = T^[ni] = T,„j r,„2 • • • T,„^ , 
where ruk = —1, or 1, this Hamiltonian reads : 



i7(4) = _ lr(2)(_i 1) 

+ -lr(3)(-i,o,i) 



(8) 



+ ij(rW(-i, 1,-1,1) - tW(-i,o,o,i) 

_ ir(4)(_;L,-l,l,i; 



This expression for the effective Hamiltonian needs to 
be written in a spin only notation. Following Ref. 0, 
one can map the singly occupied subspace of states of 
the Hubbard model onto the states of a Hilbert space of 
interacting 5' = 1/2 spins. The correspondence is: 



Hubbard space 



Spin — space 



site i 



(9) 



site i 



The expression of the spin Hamiltonian ui^^ , acting on 
the spin space, is derived from the Hamiltonian acting on 
the occupation number subspace, and is: 



(6) ^ 



1 



N 



3 

2^ ^ 

ai,Q2,'-' -a^v— \/ — 1 



where a^^} is the af Pauli matrix associated with site p. 

A more detailled discussion of this mapping is given in 
Appendix^ To order t{t/U)^, and dropping constant 
terms, we recover the results found using this method in 
Ref. I3| and found earlier, via another methodi: 

<i,j> <i:'h> <i,j3> 

<i,j,k,l> 

-(s,-Sfe)(s, • SO], (11) 

where j, j2 and J3 are respectively first, second and third 
nearest-neighbors of i, and {i,j,k,l) denotes the four 
spins that form an elementary square plaquette circu- 
lating in a clockwise direction. The coupling constants, 
homogeneous with an energy t, are expanded to 3'''^ or- 
der polynomials in t/U, giving Ji = it'^/U — 24i'*/f7^, 
J2 ^ J3 ^ U^/U^, and Jc = mt^/U^, as in Ref. 0. 
Since in what follows we only consider Hi^^ to order 
fc = 4, we henceforth use Hg as a shorthand notation 
fori/f'. 



(10) 



B. Staggered Magnetization Operator 

The Hubbard ground state wave vector, |0)h, ex- 
pressed in the effective theory, e*'^|0)H = |0)s, has a 
unique value in the SO subspace. However, it is im- 
portant to note that |0)s is not simply a projection of 
|0)h onto that spaced. In performing the transformation 
the particle excursions perpendicular to the SO space are 
taken into account in the effective theory by the non-local 
exchange integrals. The vector |0)h is therefore rotated 
by e'"^ to lie entirely in the SO subspace. Similarly, phys- 
ical quantities in the effective theory are not the expec- 
tation values for operators calculated with the projection 
of the vectors into the subspace. Since |0)h = e~*'^|0)s, 
the expectation value of an operator Oh in the origi- 
nal Hubbard model can be computed in the state |0)s 
as long as the transformed operator Os = e^'^One"''^ is 
used^i^i^°i"i^^i"i^^. In other words. 



(0) = 



h(0|Oh|0)h ^ s(0|Os|0)s 
h(0|0)h s(0|0). 



(12) 



These operators Os may differ from the expected form 
in a phenomenological magnetic model constructed 
uniquely in the SO Hilbert space. We focus here on the 
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operator for the staggered magnetization (magnetic mo- 
ment) for the Hubbard model, M^. We show that, when 
considered in the effective theory, the magnetic moment 
is not the same as the Heisenberg magnetic moment oper- 
ator M^. We henceforth use the tilde symbol to annotate 
what an operator, Os, would be in a SO model with no 
relation to an underlying Hubbard model. We define the 
conventional staggered magnetic moment operator, M^, 
that lives in the unrestricted Hilbert space of the Hub- 
bard model as Mj^ = (l/^)E^(^^,T " n.^Oi-'^Y- We 
consider a square lattice of size Lx x Ly = N, with sites 
labelled i£[l,N]. 

The unitary transformation on ( = 

e^'^Mj'je"*'^ ) is performed using the Hausdorff formula, 
as in Eq. ©, leading to the expression for a new opera- 
tor in the SO spin subspace. Such a calculation can 
be achieved with the commutation relations between 
and T_i, Tg and Ti: 



N 



hi. 



(13) 



(14) 



(15) 



where we henceforth use the notation Tm to emphasize 
the similarities between the expression for these new op- 
erators and the original kinetic operators T^, and to 
point out that both and Tm increase the number of 
doubly occupied sites by m. 

The spin Hamiltonian in (|11|) is the exact expansion of 
the original microscopic model in obtained by consid- 
ering terms in S^''^ ( Eq.(7) ), up to 3'"'' order in [t/U). 
It is important to realize that an odd power in t cannot 
appear in the spin expression of i?!'^' ( Eq.(6) ) for any k 
because we are at half-filling and all states are singly oc- 
cupied in the low-energy theory. In the case of staggered 
magnetization, transforming Af^ and retaining terms in 
up to third order in {t/U) would lead to a third order 
power in t/U in the expression of Mj'. However, the third 
order contribution in t/U to MJ' evaluated within the |0)s 
SO ground state would vanish since we are at half-filling. 
Hence, we only need to keep terms in S'^^^ up to second 
order in [t/U) (i.e. the terms proportional to 1/C/ and 
1/C/^). Doing so, we obtain the following expression for 
the effective staggered magnetization operator in the SO 
subspace: 



Ml 



u 



. 

The linear term only contributes when the expectation of 
higher powers of are computed in the |0)s SO ground 
state. Then, linear terms can combine to give an overall 
{t/UY contribution. 

Using the same procedure as above in obtaining 
Eq. ifTT)) (see Appendix ^ and restricting ourselves to 
expectation values of the first power of magnetization in 
the SO subspace, we find the expression for the magne- 
tization operator in terms of spin 1/2 operators: 



\ i <%3> J 

(17) 

This expression for contains a correction compared to 
the standard staggered moment operator in a Heisen- 
berg model. 



(18) 



This is a consequence of the fact that the original Hub- 
bard model contains electron mobility, or charge fluctu- 
ations, where particles are allowed to visit doubly occu- 
pied sites. The magnetic moment of the ground state has 
therfore non-zero contributions coming from high-energy 
configurations with doubly occupied sites. Within the 
large— ?7 limit, hopping is highly correlated and limited 
to sequences taking the system between two configura- 
tions in the SO subspace^. When represented in the 
effective theory this particle mobility gives rise to addi- 
tional quantum fluctuations over and above the quantum 
spin fluctuations of the S = \/2 spins around a Neel or- 
dered state. Hence, in calculating the magnetic moment 
in the effective theory one must use the operator and 
not Mg^, the latter being used in phenomenological stud- 
ies dissociated from a parent high-energy Hubbard-like 
fermionic niode&i^iiSiiLiSiia. 



C. Alternative Derivation 

As an alternative and possibly more physically trans- 
parent method to obtain the transformation of the stag- 
gered magnetization Mj' above, and to help in the phys- 
ical interpretation of the result, we add a conjugate field 
/ijj to the Hubbard staggered moment. 



Hi,^ Hn- hi ^(n,,T - ^a)(-l)' 



(19) 
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and repeat the unitary transformation calculation start- 
ing back at Eq. 0. We find 



0+2 



(20) 



<»,i> 



which satisfies the following relationships 
hm - — — # = MJ and lim 



_o N dhl 



with Mgt given by Eq. (|17|l and where 



(21) 



hi ^hl{l- 22— 



Aft= (l-2z— )Mt, (22) 



with z the nearest-neighbor coordination number. This 
result gives an alternate point of view (and distinction) 
of the above relationship between the SO, AI^, and Hub- 
bard, Mg, magnetic moments. The moment is the 
response to an effective microscopic conjugate field, hi, 
that is renormalized (reduced) compared with the "ap- 
plied" staggered field. This renormalization of the 
staggered field offers another interpretation for the ad- 
ditional "amplitude" fluctuations arising from the finite 
electron mobility. Since the weight of doubly occupied 
states become more important in the effective theory as 
t/U is increased, the local microscopic "spin holding" 
staggered mean-field as well as Mj{t/U) decrease with 
increasing t/U. As discussed further below, this correc- 
tion corresponds, to order (t/U)'^, to the finite t/U re- 
duction of the spin-density wave amplitude found in the 
Hartree-Fock solution of the Hubbard modeP". 



III. RESULTS AND CONSEQUENCES 

In this section we test the accuacy of the transforma- 
tion from to Mj' through exact diagonalisation of 
small clusters, the behavior of Mj' and Mj' are com- 
pared, as a function of t/U in the thermodynamic limit, 
using spin wave calculations. 

A. Behavior of Small Clusters 

As there is no broken symmetry for small systems, we 
calculate and its SO counterparts, Mj ., and M2 ., 

defined by 



Mi 



{{Ml) 



t\2\ 



and M. 



{{Ml) 



(23) 



where a G {H,s}. (Note that the canonically trans- 
formed (M^)'^ is not the square of Mt in Eq. itT^ l. 




0.0001 



FIG. 1; Difference between e''^|0)H and |0)s restricted to the 
singly occupied states. The result is compared to a {t/U)'*' 
line obtained by fitting 5 in the range t/U G [0,0.05]. The 
calculation is done on a 2 x 3 lattice with open boundary 
conditions in both the x and y directions. 

For small lattices, of size x Ly, the ground state |0)h 
and |0)s of i/n a-nd Hg can, respectively, be determined 
exactly. We find by direct inspection that the unitary 
transformation, e'"^, applied on |0)h, indeed decreases the 
spectral weight of configurations with doubly occupied 
states. As an overall measure of the quantitative agree- 
ment between e*'^|0)H and |0)s and of the accuracy with 
which the doubly occupied states are eliminated from 
|0)h, we plot in Fig. 1 



(5 = ^|(n|e''5|0)H-(n|0), 



(24) 



where the sum is carried over all 2(^--^«) singly occupied 
states. Here a system of size 2x3 with open boundary 
conditions was considered. The overlap between the two 
state vectors diminishes as t/U increases, with a differ- 
ence that is roughly proportional to (t/U)*, the order of 
the first terms neglected in the calculation. 

In Fig. 2 we show results for Afjg, M^^^ and M2S 
for a 2 X 4 system. The full curve (circles) shows re- 
sults for exact diagonalization of the Hubbard model, 
Afjjj, which should be considered as the reference data. 
One can see that M^^ is a decreasing function of t/U 
at small t/U, as expected on physical grounds and as 
found in previous exact diagonalizations'^^ and series ex- 
pansionSi. The dot-dashed curve (rhombuses) shows the 
result for Af|j,. While there is a quantitative difference 
between the two results, one finds that the two sets of 
data share the same slope, as {t/U)"^ and that their 
difference (not shown) scales as {t/U)'^ for small t/U. 
The dashed curve (squares) shows the t/U dependence 
of the magnetic moment calculated from M^^ and |0)s. 
Contrary to the exact result for M\ ^ and the SO result 
Afj g, M2 s increases with (small) t/U, and never has the 
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FIG. 2: t/U dependence of the staggered magnetization Mj g, 
jj and M2,s for a 2 x 4 lattice. 

correct limiting small t/U behavior. Simply calculating 
the staggered magnetic moment, as defined in a Heisen- 
berg model, is therefore qualitatively incorrect when the 
low-energy Hamiltonian includes higher order corrections 
in t/U . On the contrary, when the correct SO operator 
is used, the result is not only qualitatively correct, 
but the difference between the exact Hubbard result and 
the SO result is less than 1% for t/U = 0.1. This sug- 
gests that {U/U)^ = .026, with U the half-bandwidth, 
gives an estimate of the error on the staggered moment 
in the SO theory. Since the mapping between the Hub- 
bard model and the effective theory is not size dependent, 
we expect this accuracy estimate to roughly apply in the 
thermodynamic limit. 



0.305 - 



0.3 - 



s(oim:io), 

-»(0|Mt|0), 



0.02 0.04 ^ 0.06 0.08 0.1 



FIG. 3: {t/U) dependence of (Mj") and (AfJ) in a Holstein- 
Primakoff calculation of Hs to order 1/S. 



Aj: = AS{Ji -J2~J3~ JcS^ + (J2 - JcS^)iyj: 



Br. 



Ho = -mSHJi-J2-J3-^JcS^) 



and 



= 2*-*^°'^^^^'' ^ cos(fcj^)) , 
ly^ = cos{kx) cos{ky) , 

Tj: = -(cos(2/cj;) -I- cos(2fcj^)) 



(27) 



(28) 



Defining 



Bl 

k 



(29) 



B. Thermodynamic Limit: Spin Wave Calculation obtain the standard resulliSi 



As, in the absence of boundary effects, and dif- 
fer only by a multiplicative factor (see Eq. 122() ) , one can 
estimate the effect of this factor in the thermodynamic 
limit within a spin wave approximation. In this case, the 
spins operators are written in terms of their bosonic exci- 
tations through a Holstein Primakoff^^-^^ 1/S expansion 
of this bi-partite Neel ordered lattice: 



Sublattice a 
Sf = S - a\ai 
S+ = V2S[a, + OiS-^)] 
Sr = V2S[al + 0{S-^)] 



Sublattice b 



SJ- = V2S[bj + 0{S-^)] 
S+ = V2S[b] + 0{S-^)] 
(25) 

Transforming in Eq. 111(1 in reciprocal space, one 
obtains to order S 



H = Hn 



where: 



T.^k 



blbr 

k 



iibi 

k k 



(26) 



(30) 



We show in Fig. 3 the results for {M^) and {M^) cal- 
culated to order 1/5* in the Holstein-Primakoff formula- 
tion of the Hamiltonian in Eq. (|ll|l . The data show 
qualitatively the same behaviour as for the exact diago- 
nalization (see Fig. 2): a positive trend at small t/U for 
the moment of the SO model and a negative trend 
for the transformed moment M^. Even though the ring 
exchange term is of order S"^ larger than the bilinear ex- 
change terms, a calculation that would keep boson oper- 
ators beyond quadratic order is apparently not required 
to get the correct qualitative trend of vs t/U . 

These results have several immediate consequences. 
We conclude that the increase of the Neel order parame- 
ter in the presence of ring-exchange ^^^iSiiLiSiiS is due to 
the use of Ml , which neglects the renormalization factor 
(1 — 2zt^ /U"^) coming from charge quantum fluctuations. 
Further, we note that this renormalization factor is, to 
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order (t/U)'^, identical to that reducing the spin-density 
wave ampUtude in a Hartree-Fock solution of the Hub- 
bard model22i. 

We note here that it should not be construed that all 
quantities measuring the strength of the magnetic cor- 
relations need to be a monotonously decreasing func- 
tion of t/U. For example, when considering the three- 
dimensional Hubbard model, where the Neel tempera- 
ture Tn is nonzero, one finds that T^/J = 0.9575 in 
the Heisenberg limit^^, where J = At'^/U is the nearest- 
neighbor Heisenberg exchange. Normalizing Tn by the 
scale t, we have T-^/t = 3.83{t/U) in the Heisenberg 
limit^LSa. Hence Tf^/t is a non-monotonous function of 
t/U, first increasing as U decreases from the Heisenberg 
limit, and then decreasing as the weak-coupling limit 
t/U ^ 1 is approached. However, this behavior comes 
the fact Tf^/t is controlled by the spin stiffness which 
scales with t/U in the opposite manner to the zero tem- 
perature order parameter, in the strong coupling limit. 
The spin stiffness is controlled by J but the magnitude 
of the antiferromagnetic order parameter is not, as can 
be seen trivially in the Heisenberg limit where it is in- 
dependent of J. More to the point is the observation 
that even for a relatively large U/t = 10, Tj^/t is already 
25% below the Neel temperature that would be predicted 
by the Heisenberg model^^-. Here, the charge fluctua- 
tions lower TnccI of the Hubbard model below that of the 
corresponding limiting Heisenberg model. In the con- 
text of the work presented here, it would seem possible 
that a numerical calculation on a three-dimensional effec- 
tive spin-only Hamiltonian to order t{t/U)^ that neglects 
charge renormalization would give a Neel temperature 
that actually increases even faster than the Heisenberg 
T^/t = 3.83{t/U) and definitely faster than Tn, obtained 
numerically for the Hubbard model. 



model and that from the spin-only theory is of order 
{4t/U)^ where the power 4 is the first power that is ne- 
glected in the derivation of the low-energy theory. We 
also showed that even though the ring exchange term 
is of order larger than the bilinear exchange terms, 
a calculation that would keep Holstein-Primakoff boson 
operators beyond quadratic order is apparently not re- 
quired to get the correct qualitative trend of AI^ vs t/U. 

Finally, note that charge fluctuations should also lead 
to amplitude renormalization factors for magnetic order 
at other wave vectors or for other order parameters such 
as dimerization. Renormalization factors for other ef- 
fective models, such as the spin model coming from the 
three band model of the Cu02 plane, and models that in- 
clude second, t' , and third, t" , nearest- neighbor hopping 
termsS^, are also open problems. 
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APPENDIX A: DERIVATION OF THE SPIN 
HAMILTONIAN 



IV. CONCLUSION 

We have shown that transforming the Hubbard model 
into an effective spin only theory leads, for t/U away 
from the t/U ~ Heisenberg limit, to a new source of 
quantum fluctuations that reduces the staggered magne- 
tization. Indeed, short-range charge fluctuations renor- 
malize the order parameter by a factor [1 — 2z{t/U)'^], 
depending on t/U, which is independent of the spin-only 
quantum fluctuations. This factor insures that increas- 
ing the charge mobility reduces the overall stability of the 
magnetic phase at t/U > 0. This is despite the decrease 
in long range zero-point spin fluctuations which, when 
taken alone, suggests that the antiferromagnetic order 
parameter should become larger as t/U increases from 
the Heisenberg limit. It would be interesting to find out 
whether this separation of charge and spin fluctuations is 
maintained to higher order in the perturbation scheme. 

As a quantitative guide for the validity of the strong- 
coupling expansion, we also checked on small clusters 
that the difference between the result from the Hubbard 



Eq. ||SJ) gives the expression for the third order expan- 
sion of the Hubbard model in terms of Tm operators. 
Eqs. © and H10|) introduce the mapping between the 
spin S" = 1/2 operator Hilbert space and the singly oc- 
cupied subspace of the Hubbard model. In this mapping 
the calculation of the trace represents a somewhat subtle 
part because of the anticommutation relations between 
the different fermionic operators. As an example we de- 
rive the complete expression of the spin Hamiltonian to 
the first non zero order. We start with: 

= _^T(2)(_l,l), (Al) 

that is: 

- ^ V (h -C^ C- 71 -'^ 

il,»2,jl J2 

® ("ii .<y-i cli .a, Cji ,0-1 hj, ^ff,^ . ( A2) 

Since we work in the singly occupied subspace {{V) — 0), 
the two electronic processes that first increase (Ti) and 
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then decrease (T_i) the number of doubly occupied sites 
have to be performed between the same 2 sites, which 
imphes for Eq. l|X2)l that: 



J2 

i2 



ji 



(A3) 



Defining the fcrmionic orbitals as: 

nil, nix, "2T>"2i, • • • , nNi;,nMi 

1=* Site 2"'* Site Af**" Site 

For example, we have for two sites: 



(A4) 



ci,t|t,T> = ci,t|l, 0,1,0) = 10,0,1,0) 

c2,TlT,Ti) = c2,t|i, 0,1,1) = -|i, 0,0,1) 
c2alT,Ti) = C2,x|i, 0,1,1) = |i, 0,1,0) 



|o,T) 

-|T,i) 
IT,T) 

(A5) 

the minus sign coming from the odd number of occupied 
fermionic orbitals occuring before the one specified by 



Ci_CT or cj^. We use now the t, i notation for simplicity. 



but it is important to notice that it represents fermionic 
orbital occupancy. The symbol |0) represents a site that 
is empty for both its orbitals t and j. It follows that: 



TillA) 



ri|T,i) = i|o,iT) + tUT,o) 



T-i|0,iT) = 
T-i|Ti,0) = 



-tU,T) + t\hl) 
-<U,T) + t\hl) ' 



(A6) 



(A7) 



and finally: 



T-iri|t,i) 

T-iTili,!) 



2tmti) - \u)} 

-2tMIT,i) -U,T)} 



We can then rewrite ljA2|l as: 

H'^'^ = -^(IT,i) -U,T))((T,il -a,TI) 



(A8) 



(A9) 



This form makes it easier to calculate the trace in (jlOf) 
and gives: 



U 



-2^ 
U 



so that: 



Tr(aiV|i/(2)) = -2- 



(AlO) 



(All) 



(A12) 



(A13) 



(A14) 



<i,3> 



<hj> 



(A15) 



recovering the well known result of the nearest neighbor 
interaction coupling constant Ji : 



it' 



(A16) 



The same method can be applied to calculate the expres- 
sion of the spin Hamiltonian up to order t{t/U)^. As in 
Ref. 0, we used a program for the general construction of 
Hi^'> for fc > 4. To order k ^ 4, this leads to Eq. ((TTJ- 

<ij> <ij2> <i,j3> 

<i,j.k,l> 

-(S, •Sfe)(Sj • SO], (A17) 
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This means that in the large U/t limit, at some finite tem- 
perature T, antiferromagnetic correlations must increase 
with increasing t/U, going from zero staggered magnetiza- 
tion to a finite value as one crosses the phase boundary. 
(See Ref. l22|). 



